Copy — 


SUBTASK REPORT (LET 





PROBLEMS AND _ POSSIBILITIES 
OF HIGH ALTITUDE 
WORLDWIDE MAPPING 


Development of Mathematical System 
Geometry Models 


Task1, Subtask C 


Itek Data Analysis Center 


Prepared by 







Approved by 





10 December 65 





Itek Corporation 


HS iriaeialbbse 
TACENTOREYAOLE 
: ONTROD=S46tei 
1 Accordance wit & O. 12958 : 
n_NOV 26 1997 3 





CONTENTS 


Status of Task. . . . .« « « © © © «© © © « 


Introduction . . . . . « © © «© «© «© © «© « 


Selection of Orbital Model. . . . . . 6 © « » 


1.1 
1.2 
1.3 
1.4 


Introduction . . . . .« « «© «6 «© © « « 


Notation and Definitions . ..... . 


Rectangular Coordinates XYZ of Body .. . 
Coordinate Systems and Transformations . . 


Photogrammetric Considerations . . 


2.1 
2.2 
2.3 


System Configuration . . . . . . 2. « » 
Projective Relationships . . . . .« .« « « 
Camera Orientation. . . . . «©. « « «© « 


Constrained Solution . . . . . +. +6 «© « « « 


3.1 
3.2 
3.3 
3.4 
3.5 


Introduction . . 2. « « «© «© © © «© © « 
Least Squares Solution. . . . «© + « © «© 
Application Bote Ne, ah eS 
Error Propagation . . . . . 2. « « -« 
Auxiliary Data . .. .. 2. 2. «© «© 


Empirical Solution . . . . 2. « « «© «© « « » 


Experimental Design for an Error Model .. . 


3.1 
Dee 
5.3 
3.4 
3.9 


Introduction ...... 

General Remarks ..... . + « « 
Covariance Methods at Ye ce i Sah ca Sr Ge 
Applications and Analysis By dy ea ces te 8 
Relative and Absolute Accuracy Considerations 


Appendix A — Construction of Fictitious Model. .. . 
Appendix B — Detailed Formulation. . ..... . 
Appendix C — Inversion of Large Matrices .... . 





Page 


i as 
eS. 


i i 
> pet 


hh do es 


i 
NN Ee BP BE OY SEP eS Oe ee 


i] 
an an 


4 ‘ ‘ 
> 
rary 


e 
s 

AWDFrP nHHMK HHH F&F WOwWHWWwWwW NNND Bee ep 
‘ 


' 
a 


SES Vrs" 
TALENT-KEYHOWE 
CONTROL: -SYSTEM.OML 


iii 





FIGURES 


1-1 Orbital Elements. . . . «© © «© «© «© © © «© © © « 
1-2 Siderial and Terrestrial Systems . . .. . . «6 « « 
2-1 System Configuration . . . «© «© « «6 «© © «© «© © « 
2-2 PanGeometry. . . .« © © © «© © © © © «© © © « 
2-3 Geocentric and Local Coordinate Systems ..... . 
5-1 Error Ellipsoid .. : 

5-2 Standard Errors ox in Position (1), Ground X (2), Y (3), and 
Z (4) Coordinates as Functions of Standard Error in 

Photo Coordinates (z).) 4 «4 2s. Se enw 


TABLE 


5 ~ 1 NS) che matic For mulation e e ° ° e e e ° . ° e . 





Page 


9-9 





STATUS OF TASK 


This report details the philosophy mathematical exercises required in the 
generation of the MSGM. The actual Error Model Program has been written as 
of the publishing of this report; however, it has not been completely debugged. 

It has been deemed necessary to have realistic data inputs for the program from 
actual operational materials; thus, the actual massaging of the program will not 
be completed until the first PG Mission has been realized. It is conceivable that 
the computer program will receive its final polishing in the actual analysis of the 


PG System output. 


Furthermore, it is intended to incorporate both the analytical philosophy and 
mathematical techniques in the actual data reduction and mapping operations. 
This effort is to be carried forth in the continuation of this subtask and in Task 1 
Subtask D, and Equipment Procedures for Calibrated Panoramic Data Handling. 








Fig. 1-1 — Orbital elements 








INTRODUCTION 


The objective of this report is to devise an appropriate error model, whereby 
the analysis of a satellite borne camera system might be facilitated. 


The camera system consists of three prime sensors, viz., a nominally vertical 
frame camera and two inclined panoramic cameras. The latter sweep cross track, 
and are configured in such a way that symmetric convergent stereo coverage is 
furnished. This coverage is common with a portion of the stereo-overlap that is 
obtained with successive overlapping frame exposures. 


The selected error model will consider four exposures; two panoramic and two 
frame, which are constrained to orbital and auxiliary data. This is sufficiently general 
to consider all combinations of frame and panoramic exposures, subject to constraints, 
with or without stereo overlap, and may be readily extended to complete blocks and 
strips of photographs. 


The effects of various errors in the parameters of this four photoblock will be 
obtained through an analytical experimental design, that parallels the statistician’s 
controlled experiments. For this design exact fictitious data is generated, and the 
effects of varying the accuracies of different parameters (individually and in com- 
bination) on the end result are obtained by variance-covariance analyses. 


It has been suggested that the reduction of the panoramic photography might be 
accomplished with greater expediency and precision by using an empirical fitting 
process. This suggestion is based on the undefined internal geometry of many 
panoramic cameras and the instability of this internal geometry. Without a precise 
knowledge of the interior orientation, an analytical solution must fail. 


Consequently this report will consist of discussions around the following topics: 


Selection of an orbital model 
Photogrammetric considerations 

A generalized least squares Solution 

An empirical fitting of the panoramic material 
The experimental design 
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Although the mathematical development is simple, some rather extensive algebraic 
expressions result. Consequently, the step by step descriptions of the fictitious model 
data computations, and the specific algebraic expressions for various terms, are 
contained as appendices to the main text. 
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t = Epoch, the time at which the body passes through the pericenter 
t = the time at which the body is at some point P 


p = Universal gravitational constant 
n = mean motion of body: n= u¥, a3” 
M = mean anomaly: M = alt - 7) 


E = eccentric anomaly, defined by M = E - e sinE 
® = celestial longitude of body, angle XZP- 
6 = celestial latitude of body, angle CSP 


1.3 RECTANGULAR COORDINATES XYZ OF BODY 


The coordinates of P with respect to the inertial system are given by 


xX cos @ cos 8 cos ® 
Y!|=r|cos8$]/=r |cos @ sin ® (1) 
Z cos y sin @ 


where a, 8, y, are the direction angles of the vector r. The direction cosines may be 
derived from the triangles PNX and PNY (Figure 1-1) as 


cos @ = cos 2 cos (w + f) - sin Q sin (w + f) cos I (2) 

cos 8 = sin 2 cos (w + f) + cos Q sin (w + f) cos I (3) 
and 

cos y = sin (w + f) sin I (4) 


Substituting (2), (3), and (4) into (1), expanding cos (w +f) and sin (w +f), and using 
the identities 


r cos f = a(lcos E - e) (5) 
and 
1/2 
rsinf=a(1-e?) sinE (6) 


reduces (1) to the form 


sin E 


X] [Ay By fe (E) - ’ 


¥|= [Ay By 
Z} LA, Bz 








Ax Bx 
in which the matrix | Ay By |, is given by 
Az B 


Ay By cos 2 - cosI sin Q eee | 0 
1/2 
Ay By|=a|sin2 cos Icos2 0 (1- 2)” 


A, B, 0 sin I 


(8) 


sinw cosw 


This may be expanded to yield 


Ax = a (cos 2 cos w - sin Q sin w cos I) (9) 
Ay = a (sin 2 cos w + cos 2 sin w cos I) (10) 
A, = a (sin I sin w) (11) 
B, = -a (1 - e?)'” (cos @ sin w + sin 2 cos w cos I) (12) 
By = -a (1 - 2)” (sin @ sin w - cos 2 cos w cos I) (13) 
Bz =a (1 - e*)*” (sin I cos w) (14) 


The positional coordinates, given by (7) are in an inertial system, with reference 
to which the orbital elements are given. In order that they may be referred to some 
other coordinate system, transformations will be necessary. 


1.4 COORDINATE SYSTEMS AND TRANSFORMATIONS 


For mapping purposes, a system of terrestrial rectangular coordinates, Xp Yr 
Zy is convenient, especially for points, such as vehicle positions, which are remote 
from the earth’s surface. 


This system has its origin at the center of gravity of the earth, or of the reference 
spheroid, and is oriented such that the +Zy axis is directed to the mean north pole, * 
and the mean meridian of Greenwich{ lies inthe X; ZT plane. This system is 
related to the geodetic longitude and latitude, A, ~’, and the geocentric longitude and 
latitude, A, v, according to 


* As defined by Service International des Latitudes. 
+ As defined by Bureau International de 1’ Heure. 








xX (N + H) cos yg’ cosa cos y cos X 
Y} =| (N+H) cos’ sind |=rjcosgsinA (15) 


T [4 -e’)N+H] sin y’ sin 9 


where N = radius of curvature in the prime vertical 
H = elevation of point above reference datum 
r = geocentric radius of point 
e = eccentricity of reference ellipsoid 


It is to be noted that this coordinate system is fixed with respect to the earth’s 
surface. 


For points which do not rotate with the earth as, for example, discrete orbital 
positions, it is desirable to define a coordinate system independent of the earth’s 
rotation. The sidereal system is such a reference framework. 


In this system, the Z, axis corresponds with the instantaneous pole, and the 
Xs axis with the apparent vernal equinox, y. The relation between the T and S 
systems is illustrated by Figure 1-2, the transformation being: 


x x 
Y} =MM,/Y] | (16) 
ZJs ZIT 
where 
cost’ -sint’ 0 
M,=|sint’ cost’ 0 (17) 
0 0 1 
and 


cos é sing sinn -siné cos7 
M,=| 0 cos7 sin 7 (18) 


sin -cos — sinn cosécos7 
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Fig. 1-2 — Siderial and terrestrial systems 
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in which t’ is the Greenwhich apparent sidereal time in UTI,* and &, kK are the 
coordinates of the instantaneous pole.t+ 


The matrix M, may be considered to be a unit matrix for practical purposes. 
The component angles , 7 are less than one second of arc. Furthermore, for short 
periods of time the changes in &, 7 are so small that Zy7 and Zs may be considered 
to be coincident. 


*UTO is observed time; UTI is UTO corrected for motion of the pole; UT 2 is 
UTI corrected for seasonal variation in the earth’s rotation. Bulletin Horaire No. 4, 
Series 4, 1955. 
tMelchoir P. 1954. Observatoire Royale de Belgique. Monographie No. 3, 
Service International des Latitudes. 
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2. PHOTOGRAMMETRIC CONSIDERATIONS 


2.1 SYSTEM CONFIGURATION 


Consider two overlapping frame photographs, nominally vertical, exposed at 
times t, and t,. Let two symmetric convengent panoramic photographs, exposed at 
times t, and t, overlap within the stereo model formed by the frame photographs. 
This is illustrated by Figure 2-1. 


2.2 PROJECTIVE RELATIONSHIPS 


For convenience, let the panoramic photography be transformed into equivalent 
frame photography, according to Figure 2-2. 


Define the principal point as that position on the cylindrical panoramic format 
intersected by the camera z-axis when that axis is at the mid-point of the scan. 
The equivalent frame photograph is tangent to the generator passing through this 
point. 


An image point x’, y’ on the pan photograph may be transformed into the equivalent 
frame coordinates x, y. Consider the case for diapositives, and define the +y~-axis 
in the flight direction, and the +z-axis upwards. This right hand system is illustrated 
in Figure 2-2. 


The scan angle ranges from -a through zero to +a@ as the optical axis scans from 
+x through zero to -x. Let the scan rate, &, be signed positive for a forward looking 
camera. Then 

(x’ - xj) = - fa 
(y’ - yo) = (y - yo) cos a (19) 
(x - x9) =-ftana 


with image motion compensation given by 


soa SS ein SC sine 
(Z-Z)) & ; 
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Pan model 


Fig. 2-1 — System configuration 
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Fig. 2-2 — Pan geometry 
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where 
ow ML (20) 


then 

(y - yo) = [ly’ - yg) + C sin a]/cos a 
It is to be noted that 

& (ty - ty) = ay = (x’ - x)/- £’ 


The aft looking camera scans in the opposite direction, so that @ is negative, and 
a@ goes from + to - values as the scan goes from -x” to +x”. As before 


(x” -Xg =-f’a 
(x -x)) =-f” tana 


(y” - yj) = (y - yp) cos a 
noting that @ is + when x” is -; and the image motion is given by 


V a: Care . 
(Z- Z,) % ‘Sina= -C sina 


where & is unsigned, or +C sin a if @ is signed. Then 
(y - y.) = lly” - y%) + C sin alos a 
provided that @ carries the appropriate sign. 
Xj - Xo 
Any image vector | y; - yo} in the equivalent frame coordinate system, is related to 
= f Xj = Xo 


the corresponding object space vector, | Yj - Yo}, according to 


Zi - Zo 


4 
ar 
5 
Et 


(21) 








from which the projective relationships 


wg eg [Pa Ki Xo) + mye (Vi - Yo) + my (Zj - Za} (5) 
At Pao ™3; (Xj = X,) + M39 (Yj - Y)) + M33 (Zi - Z) 


and 


) ) 


nee mas (Kj = No) + May Yh = Vo) ¢ may Zi - Zo} (23) 
i 0 M3; Xj - Xp + M39 Y; - Yo + M33 Zi - Zo 
are obtained. 


In these expressions f is the camera focal length; Mjk is an element of the 
orientation matrix Mj; X 9, Yo, Zp are the coordinates of the exposure Station in the 
XYZ system; Xj, Yj, Zi, are the coordinates of the ground point corresponding to the 
image point xj, yi, f. 


2.3 CAMERA ORIENTATION 


In conjunction with the panoramic cameras, horizon sensors are utilized, which 
yield initial and final values of the camera roll and pitch with respect to the local 
horizon. At these specific times, the local horizon is equivalent to the tangent plane 
at the nadir point, which is defined as the point at which the position vector F to the 
air station intersects the reference ellipsoid. The longitude, A», and geocentric 
latitude, $o, of the nadir, as illustrated by Figure 2-3, are determined from 


Xo cos $y COS Ay 
Yo =r {cos By sin Ag (24) 
Zs Sin Bp 


according to 


"1/2 
Ay = cos? [(xox} ee a4 | for Xo = Yo (25) 
or 
int [(v,)ox3 + x3) | < 
Aq = Sin [ Yo\Xo + Yo do for Xo <= Yo (26) 
and 
x -1/2 
by = tan” (2x3 + ¥}) | (27) 
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Fig. 2-3 — Geocentric and local coordinate systems 
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With reference to Equation (21), the orientation matrix rotating the ground coordinates 
into the photo system may be expressed as 


cos Kj sin Kj Ol! cos j Q -sin Oj 1 0 0 
Mj = }-sin Kj; cos kj Oj) 0 1 0 0 coswj sin 4; (28) 
0 0 1j{sin ?j 0 cos gj Q -sin Wy COS Wj 


The orientation between the local coordinate system and the geocentric terrestrial 
system is determined directly from Figure 2-3 as 


cos Ay -sin Ay O|]f sin $, 0 cos ®o] f-cos Aq -sinArA, 0 
Oj= {sin Ay cos Ay 0 0 1 0 Sin A» -cOSAy 0 (29) 
0 0 1{|-cos $, 0 sin ®, 0 0 1 


Equation (21) may now be rewritten as 


Xj = MjOj [xj - Xo, (30) 
or 
Xj 7 Xp Xj ~ Xp 
yi - Yo! = [XG ]lojllwj Aol So][%0]} Yi - Yo (31) 
-f Zi - Zo 


The values of all terms in (31) pertain to a discrete time. For any specific time to, 
the corresponding terms may be evaluated and held constant. In this manner, the 
dynamic local system may be transformed into a static one, with reference to which 
subsequent computations may be performed, for any time tj, provided that appropriate 
variations of Xo, due to orbital motion, of Wj, Qj) Kj, due to capsule tumbling and of 
X; due to the earth’s rotation, are applied. 


This is accomplished through (16), computing in the sidereal system, and 
appropriately expressing variations in the angular orientations 


Qj = Wj, Yj OF Kj as ay = (arg) + Cary ty; + (a, ); ead (32) 


so that the orientation matrices la;] in (31) become 


[oj] = [oo], or ty] [or tH], | (33) 
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although information may exist which indicates thatthe aj have a specific nature, 
say periodic. In this event the appropriate expressions for the ay must be used. 


Denoting the product M,M, in (16) by Tj, Equation (30) may be rewritten as 
Xj = MjO; 'T}Xip - Koy, (34) 
or 


This is in a form that is identical with (21), so that the projective Equations (22) and 
(23) may be used to define the functions 





Gyj = Xj - Xo 

Goj = ¥j - Yo 
These two equations may be expressed in termsof the constituent variables w;, 
Djs Kjy Any ee ++ Zo. It is noted, however that some of these are functionally 


related to the orbital parameters, and to each other. The constrained solution 
must utilize independent parameters, as outlined in the subsequent section. 
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3. CONSTRAINED SOLUTION 


3.1 INTRODUCTION 


The functions given by Equation (35) express the relationship between a series of 
data consisting of observed variables, and various parameters, known and unknown. 
Those which are known, are derived from observational data and are subject to the 
errors in these fundamental data. 


On substituting the various parameters and data into (35), the value Fyj will be 
zero if, and only if all entries are exact. The problem is to determine the most 
probable values for all parameters and variables. This is accomplished through a 
constrained least squares solution. 
3.2 LEAST SQUARES SOLUTION 


Equation (35) is considered to be a function of the observed variables 
iy Yi XL, Vis XY, Yi 


where the primed and double primed coordinates refer to the two pan cameras, and 
a function of the parameters 


Xo, Yos Xo, Yo Xo, Yo f, ag e 
Wj, 9}; Kj» Q, w, Ie, n, 7; Xi Yi» Zj 
Consequently, rewrite (35) as 
Fri = fx (xy, yi, Xi, Yi, X1, YL3 Xo» Yoo £ Xo» Yoo £’, X0y Yoo £73 


Wj; Qj; Kjp 2, w, Ie, n, 7; Xj, Yip Z;) (36) 


Assuming that approximate values of these parameters are available, designated by 
a superscript °, an approximate value of Equation (36) is obtained as 


° 
? 


Fi = fy [xf yf (xp, (yp, |... XP YP 27 (37) 
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The true values of these items is obtained by applying a correction to these approximate 
values, according to 


Sh eee. 
Xj = Xj + Pk xj 


Yi = yi + Pk m3 


= (x’)° + De 3 oxi 
ae? + pk dy} 
xj = (x{)° + py dx] 
yii= (yf) + pk byf 


i Ae 
Yo= Yor Qk SYo 
f= f°+q) of 
x9 = (xg) + qh xy 
yo = (yo)? + ak Syo 
f’ = (£1)° + gf of’ 
= (x7)°+ qh dx 
ee cae: 
f=")? + ah Of” 
wes wp + ay Su; 


Qj = OF + at 59; 
Kj = Kj + qi}? 57; 


Q = 2° + qi 62 


w= w+ git bw 
I= [°+qy I 
e=e°+q) be 


Y¥,=Yj+ay oY; 
Zi = Z]+ qi! 6Z; 








the values 6xj ..... 652Z; being corrections, the values p,q being equal to unity or 
zero, depending on whether the element is unknown, or known exactly. 


The reduced condition equations may be expressed in the form 


where Vx = a vector of unknown residuals 
A, = a vector of unknown parameter corrections 
4, = a vector of unknown parameter corrections 
Ey = a vector of random variables 


and 


power NLS ee ere (39) 


* ~ 8 (unknown variables) 2 (xj, yi j 


B, 7 ) (F,, F.)ij 7 ) (F,, F2)ij (40) 


~ 8 (unknown parameters) 2 (all or none or x), ... 7) 


— a (Fy, Fy)ij - a(Fy, Fo)ij a 
2 a (parameters to be constrained) 2 (remaining parameters) : 


Assume that statistical estimates of say q parameters are known, and that we 
wish to constrain the adjustment to fit these estimates. Designate the statistical 
estimates of the parameters f by 


B° =| Be (42) 


having an associated covariance matrix 9g: 


078, OB,B. ..+- OByBe 
OBB, OB. ... OBBq 
9B = (43) 
OBB, OB _B2 eee 7By 
3-3 
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where 8° =a current corrected value of 6° 


Va = an unknown residual vector 
This equation is reformed as 


Wyte se 


whence Equation (38) may be rewritten as 


AV +BA+E=0 


where 


Vx 
7 = 

Vp 

ay 
Ae 

4, 

Ex 
E= 

G 

Ax, 0 
A= 

0 iI 





(44) 


(45) 


(46) 


(47) 


(48) 


(50) 





and 


B, B, 


-Iqq 9 
The conditioned solution of Equation (46) for V and 4 which minimizes 
Ss=vio'v (52) 


is required, in which 


G= (53) 
0 og 


where 0x = the covariance rnatrix of the observed variables 
The required solution is obtained from 
S=VT o-ty - aT (AV + BA + E) 7 (54) 


in which A is a vector of Lagrangian multipliers, setting 


= =0 and ve 0 
to yield 

V=ocATa (55) 
and 

-2BT A =0 (56) 


which are combined with (46) to yield 


n= -(A oAT) (BA +E) (57) 


Oe 
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and 
-1 
A=-[BT(AcAT)-'B] BT (AcA!)"E (58) 
This equation is evaluated by partitioning, the partitions being given by the following 
identities: * 


Ax Ox ys 0 


AcA! = (59) 


-1 
AcAl = (60) 
0 oR v1 
Put 
(Ay Ox AS): = Wx 
B? W,B,+og-! BP W, By 
BT(Ac Al)! B = (61) 
T T 
Bi Wx B, By; Wx B, 
and 
T - 
B, Wx Ex - 9g7'G 
BT(Ag AT) B= (62) 


B? Wx Ex 
3.3 APPLICATION 


In order to clarify the previous treatment, anillustrative example of applying these 
formulas to our problem is given. 


*The utilization of these partitions is described in Appendix C. 
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For each photograph, j, one obtains for each point i 


aF,/ox;, 9F,/ayj;, OF ,/axi, OF ,/ayj, oF ,/ax{, 9F,/ay; 
aF,/9xi, 9F)/ ayi, OF,/2xj, OF,/ayj, oF,/axj, oF,/ayf 


or 
Ay = [Ax, Ax, Ax] i (64) 
For the frame camera Axj, Axj, are zero; for the forward looking pan camera, 


Axj, Axj are zero, and for the aft looking camera Axj, Axj are zero. 


To each ground point that is known, a variance-covariance matrix may be 
assigned. For unknown ground points, coordinate values are estimated, together 
with large variances for these estimates. 


Orbital parameters are estimated from the given ephemeris, and values of the 
rotations wj, yj, Kj are determined from horizon sensors. Values for Xo, Yo, f; 
Xo) Yo {’ are obtained from calibration data. 


The end result is that all parameters are to be constrained~some extremely 
loosely—dictating that B, = 0. 


Consequently, we may write (see Equation (65) on following page) 


The various covariance matrices are: 
Ox, for the measured frame photograph images 
Ox} for the measured forward pan photograph images 
Ox i for the measured rear pan photograph images 
oXj for the ground coordinates 
oc for the frame camera constants 
oc’ for the forward pan camera constants 
oc” for the aft pan camera constants 


oa for the camera tilt angles 
and 


0, for the orbital elements 
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OF ,/2x, 
aF,/ ayy 
aF,/ df 

aF,/axj 
aF,/ ay4 
aF,/ af’ 
aF,/axf 
oF,/ ay? 
aF,/ at” 
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aF,/ 39; 
aF,/ aK, 
aF,/ aa 
aF,/ aw 
aF,/al 
aF,/ de 
aF,/an 
aF,/ar 
aF,/ 8X} 
aF,/ 8Y; 
aF,/aZj 


aF,/ aXq T 
aF,/ ayo 
aF,/ af 


aF,/ ax9 
aF,/ 2y6 
aF,/ af’ 
aF,/ ax 
oF»! ay6 
aF,/ af” 
BF ,/ 8w; 
aF,/ 29; oe 
aF,/ ak j 
aF,/ aQ 
aF,/ aw 
aF,/ al 
aF,/ de 
aF,/ an 
aF,/ aT 
aF,/ aX; 
aF,/ aY; 
aF,/ aZ; 
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It is to be noted that 0, may not be explicitly given, but that for a sequence of n 
orbital positions a 3n <x 3n variance-covariance matrix of these positions is given. 
The method of determing the covariance matrix of the orbital parameters is outlined 
below. 


Consider the parameters p = (Q, w, I, e, n, 7). By definition, the 6 x 6 covariance 
matrix of these parameters, C (p), is the expected value 


an] [an]? 
dw] | dw 
di di 
C(p) = E es (66) 
dn} | dn 
di d7 
which may be rewritten as 
E [dQ] E [da] |T 
E [dw] E [dw] 
c® = E [di] E [di] (67) 


E [de] E [de] 
E [dn] E [dn] 
E [dr] E [dz] 


These parameters are functionally related to the variables Xj, Y;, Zj, and thus each of 
the individual expectancies may be formulated as 


E {pt = a E fax} + <P & {ayy}+ <P & faz} (g8) 


i i 


By successive substitution, and factoring it is found that Cp can be expressed as 
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aQ/ ax, a2/a¥; 32/22; 
aw/ aX} aw/3Yy aw/ aZ; 


an/ ax; 20/ ay; 
aw/ aX, aw/ dYj 
ai/ ax; 2i/ aY; 
ae/aX; ae/ aYj 


ai/aX, 2i/eY, 2i/2Z; 
ae/aX, ae/a¥y de/2Zj 
an/axy an/aY, an/aZj 


at/ aX; at/ ay; at/aZ: 


C(p); = 


an/ aX; an/ ay; 
a7/ aX; 3a7/ ai 





Noting that the central factor is the covariance matrix of XiYiZi, denoting the 6 x 3 
Jacobian as J(p; X), then we may express C(p) in terms of C(Xj) according to 


C(p); = J(p; Xj) C(K,) ITB; X,) (70) 


It is to be noted that C(p)j is a variable, depending on the values of Xj and the values 
of C(Xj), which however becomes a constant if C(X;) varies in the appropriate fashion. 
This implies a dependence between successive values of Xi, as substantiated by fact. 
Thus the covariance matrix C(p) determined from m successive positions may be 
written as 


C(p) = I; X) CH) ITP; XP i=1l,...m 


6x6 6x 3m 3m 3m 3m~x 6 (71) 


which is readily soluble provided that the 3m x 3m covariance matrix C(Xj) is 
known. It is noted that Equations (71) and (68) are identical in meaning, since the 
right hand side of (71) is a summation over m points, which is in accordance with 
the definition of an expected value. 


The two condition equations for each ground point image are: 


Aj Vi+ BY 4+ Ej =0 






where 
Axi 0 
(2x6) (2x 21) 

A, = 6 I (72) 

(21 x 6) (21 x 21) 
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Now, 


Vi = [Vxi» Vy} Vxip eee ee Vr, VXir VYi; vat (73) 
By 

er a (74) 
-I 

(21 x 21) 
Ai = [oxy Syq df, 8%, Sys, Of... -- OZi]T (75) 
and 
Ej = [Ex;, Ey; Exis Ey} eevee Ez;] T (76) 


The associated weight matrix is 

oo = Diagonal [ox;, oxi, Ix{, Tes Te, FE, %, Ta, TX] (77) 
On forming the normal equations, rewritten as NA = -C, one obtains 

N = BY (Aya?) Bi (78) 
and 

C = BP (AycAT)* Ej (79) 
from which the solution for A is obtained as 

A= -N'C (80) 
The formal solution of (80) indicates that it will be necessary to invert a 
(16 + 35 + 3;) x (16 + 35 + 3;) matrix as a minimum. If the angles Wj, Pj» Kj are 
expressed in the form a = a) + a,t + a,t? +amt™, the order of the matrix is increased 
by 3mj- This is so formidable a task as to yield it impractical, even for existing 


computers. Various techniques have been developed for the inversion of large 
matrices. These are described in Appendix C. 


However, if one reconsiders the normal equations, it is found that they consist of 
two parts—one pertaining to the orbital and camera parameter data, and one relating 
to the ground point data. 
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Following the notation and derivation of Brown, the observation Equation (46) may 
be rewritten in the form 





AV + Bd + BO+ E=0 (81) 
where B, 5 refer to orbital and camera parameters 
B, 6 refer to the ground point data 
Equation (60) may be rewritten in the form 
Wx 0 07 [lAxo,AT)" 0 0 
w=(AcAT)'=|0 WwW oOl=] 0 (3)-! 9 (82) 
0 oOo WwW 0 0 6)" 
where o = Diagonal [a¢, Oe’, Oe", GO, oa;] 
O = Diagonal [ox,, OX +e ee ox;] 
Equation (82) is equivalent to 
A 0 O]fox,7! 0 OTfAT 0 O 
w={0 I ol] 0 WwW Oj/o ro 
0 0 IJLO0 owlo o1 
Similarly (61) may be rewritten as 
Wx 0 O|/B B 
N=BlwB=|BT -1 o!]0 w o|/-1 0 
BT o -ri1Lo o wjlo 1 
(BIWXB + W) (BTW,B) N+W N 
= (83) 


(BTW,B) (BTw,B+W)] LN-T New 








Similarly, Equation (62) becomes 
BF -I O][Wy 0 OJ[E 
C = BIWE = o WwW O||E 
BT 0 -I}Lo oO WILE 
or 


(BTWxE, - WE)] [C - WE 





Cz 2 (84) 
(BTWxEx - WE)| LC -WE 
The normal equations, NA = -C, are thus 
N+wn [[5 C - WE 
=- (85) 
NT NiwWié C -WE 
which for computational purposes is partitioned: 
Now! NN | OM 
T ee ee 
Be i, + W, ~ (86) 
NT | N, + W, 
Ni | Ni+W 
Let N~ =M, according to 
N+w WN (7? IM M™ 
= (87) 
BT ON.) LT i 
Since NM =I, then, noting that M = -M N (N + W)7!, one obtains 
M = [(N + W) - NN + WENT] (88) 
and 
M = [(N+ W)7!+ (N+ W)-! NIM NON + W)74] (89) 





FAR erate 









In order to determine the corrections 4, 8, it is not necessary to solve Equation (89). 
Consider 


6 M M C-WE 
J = (90) 


from which 


-6=M(C - WE) + M(C - WE) 


and 

-8 = MT (C - WE) +M (C - WE) 
Putting 

Q=(%+ Wt NT (91) 
then 

M = -MQ 
and 

-5=M(C - WE -Q(C - WE) (92) 
Since 

- [NT6 +(N+W) 6]=C -WE 
then 


- (8) = (N= W)"1(C - WE) - Q6 (93) 


It is to be noted that N + W consists of i diagonally arranged 3 x 3 matrices, and 
consequently presents no difficulties. M requires the inversion of a (16 + 3ynj) matrix. 
For the purpose of this project, this is not considered to be too large; however in 

the event of large m, j, N + W may be suitably partitioned. 








3.4 ERROR PROPAGATION 


In the final iteration 


Vx =E,-Bb-Bo 


and 


with § and 6 tending to zero, i.e., 


x Ex 
v-t 
v-# 


Now 


S=Viw, V+Viwv+vIwy 
x 


(94) 


Equation (94) which when divided by the degrees of freedom gives the unit variance 9). 


M is the covariance matrix of the adjusted camera and orbital parameters, and 


M is that of the ground po 
the covariance matrix of 


ny point g is 


Mg a (Ne + We)? + Qs Mg Qz 


3.5 AUXILIARY DATA 


ts. It is to be noted that M = (N+ W)-!+QM QI, so that 


(95) 


Suppose auxiliary data, independent of the camera system, has been collected, 


which may be expressed a 


s a function of the various parameters. As an example, a 


radar altimeter will indicate the value of nadiral distance, which may be expressed 
in terms of the orbital parameters. 


Denote the vector of auxiliary data as A, which may be written as 


Ap = Fy (Q, w,. 
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then, as before 


An = Af + VAy (97) 
The various parameters Qp = MOPS Gs iene fe ca T, ) may be expressed as 

Ap = a + bay (98) 
so that 

Van ~ Ah, 5% +--+ Ahp 5p = EAh (99) 
where 

Ean = - Fh (Qh, Og, --5-- ap) (100) 


Consequently, may write an additional series of equations: 
Va -AG=Eq (101) 


which, together with the covariance matrix 0,, may be incorporated into the previous 
solution, according to 


Vx B B Ex 
VA -A 01/5] |E, 
+ = (102) 
V -r o| {6} |E 
V 0 I E 
to yield 
(BTW,B+W+ATw, A) BIWB N+Ww+ATw, AT N 
N = = 
BTW,B (BTW,B+W)| [NT N+W 
(103) 
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and 
C -WE -AWAEA 
Cc - WE 


Similarly, auxiliary data pertaining to ground point data, (X; Y; Z;) might be used to 
exploit the relationship between auxiliary data and the various parameters. 
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4. EMPIRICAL SOLUTION 


Let it be supposed that it is desired to express the ground coordinates of points 
as a polynomial, in terms of the panoramic image coordinates, (x; - xq), (yj - yo). In 
order to do this with some discrimination, let the projective relationship be expanded. 
Consider the form 


Xj Xp -f sin aj 
Yi] = [Yo] + KiM |(y; - yo) +C sin aj (105) 
Zi Zo -f cos aj 


in which Kj is a variable scale factor, and M is the orientation matrix. 


Now, Kj is the ratio of the length of the object space vector to that of the image 
space vector, which for a tilted photograph is obtained as 


K, = H/t cos t cos ay (106) 
which may be rewritten as 


me | 
“fecost 





K [1 + a2/2! + 5Sa4/4! + 6108/6! +...] (107) 


in which the suffix i has been supressed. 


Denoting the elements of the matrix M as mjj, noting that the suffix i is not the 
same as was previously used, the component Vj, i = 1,2,3, of the last term of (105) 
may be rewritten as 


mj, [-f (a@°/3! + o°/5!-...)] 
Vy ~ K{ mj [(yi - yo) + C (a - 3/3! + a°/5! -...)] (108) 
mj3 [-f (1 - @?/2! + of /4! - of /61 +... )] 





—Fop- 92073 ur 
oe 


Substituting (107) and (108), and putting k = H/f cos t yields 

(mj, C - my f)(1 + a?/2! + Sa‘/4! + 610°/6! +...) 
(a - a°/3! + a°/5! -...) 

+mjp (yj - yo(1 + a2/2! + 5at/4! + 610/61 +...) 


Vizk (109) 


-mjg f 
This yields the odd ordered polynomial of the form 


Vi =A+B (yj - yo) +C (xj - Xo) + D (yj - yolxy - x9)? + E (xj - xo)? +... 
(110) 


Since aj = (x; - Xo). Provided that aj < 90 degrees, and that sufficient control points 
are available for a solution, an empirical fitting using Equation (110) will furnish 
ground coordinate values. 
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5. EXPERIMENTAL DESIGN FOR AN ERROR MODEL 


5.1 INTRODUCTION 


The objective of prescribing and adhering to a designed experiment is to determine 
the effects of various factors on the output data. In practice, the use of real observa- 
tions dictates that the experimental data are acquired under specified conditions. 

These conditions are then individually and sequentially varied. These classified data 
are then subjected to statistical techniques through which the effects and interaction 
of the varied conditions on the end product may be determined. 


Without real data, fictitous data may be simulated from a mathematical model, 
and treated in the same way. Although this, and Monte Carlo techniques, are frequently 
the only solution to a specific problem, the present error model lends itself to a more 
elegant and ecomonic method of analysis. This is discussed in the following sections, 
in which the following assumptions are made. 


1. That there will be four classes of ground control data, the classification being 
dependent on the precision of these data 


id 


. That sufficient ground control data is recognized on frame photography for 
the location of the exposure stations 


3. That panoramic records may contain zero or more absolute ground control 
points 


4. That panoramic records may be subdivided for mapping purposes into three 
or four segments, each containing a minimum of nine well distributed photo- 
grammetrically determined control points 


5. That the data reduction procedure consists of one of the following combinations: 


a. One frame and one pan photograph 

b. One frame and two stereo pan photographs 

c. One pan photograph and two stereo pan photographs 

d. Two stereo frame and two stereo pan photographs, in which the pan photog- 
raphy always overlaps the frame (stereo) imagery 


6. That a sequence of one or more combinations of the preceding variations may 
be adjusted simultaneously 









Sanh ot 


7. That various constraints in addition to parameters describing the vehicle path | 
may be employed in obtaining a solution 


5.2 GENERAL REMARKS 


Although the techniques of multi-factorial analysis are well known, and frequently 
used, the usual approach when using simulated data is to apply perturbations to these 
data and utilize a Monte Carlo technique. Little, or no use is made of sets of exact 
data with variations in the weighting functions for these data. 


This may be a consequence of a result being presented in the form of a variance- 
covariance matrix, rather than in the form of the easily understood residual errors. 
Provided that one has a set of exact data, the error term in the reduced condition 
equations is always zero. Consequently the vector of parameter corrections and of 
residual errors must always be zero. Consequently, no matter what weights one 
assigns to the exact input data, the resulting parameter determination will always 
furnish the same exact values. They are however, different in the sense that each 
result is associated with factors applied to the original input material. 


Varying the weights of a set of data is far simpler and more economic of computer 
time than generation random numbers to perturb data prior to a solution. Furthermore, 
the analysis of the resulting computations is simplified. 


The suggested method is to select a set of exact data, from which one may form 
the weighted normal equations using zero variances for all but one parameter or 
variable. By assigning a sequence of variances to this parameter, the independent 
effects of these precisions on the output data is indicated by the resulting variance - 
covariance matrices. By sequentially varying the combinations of weights, different 
covariance matrices are obtained, from which a complete analysis of the system can 
be made. 


In order that the various techniques might be better understood, and more readily 
interpreted, the following section is devoted to some comments on the character of 
the variance-covariance matrix. For those familiar with covariance techniques, this 
section may be ignored. 


5.3 COVARIANCE METHODS 


The covariance of two variables x;, and xj is defined as the expected value 
E {(x; - uillxj - uy)} =E {xxj)} -EkjE {xj} = 0495 pij 


where “i,j designates the mean value of Xije 
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Similarly, the variance of x; is defined as the expected value 
E {(x; . u4)* =E {xi} -(E {x;}) = of 
Assuming that the observational errors ej, of a m-variate distribution have a joint 
distribution with zero means, then oj = E {ej} and 949; pij = E {ej ej}. These are the 


general terms of ann matrix, M, which is termed the variance-covariance matrix. 


Mi, My .-- mpl [Efe} Efe,e}...E fee 


Mo, Mog «2 e Mon E {e, e,} E {e3} FE {€» ent 


mn, Mn, .-- Mpn| |E {ene} E {ene}... E {ed} 


This matrix may be normalized by transforming it into the correlation matrix, 
Pp, in which the elements Ppij are obtained as 


Covariance (xj x;) mij 
-= = 
Pij [Variance (x;) Variance (x)? 7 [my - m3)" 2 


Obviously if i = j, pyj = 1. 


For a multi-variate case, these correlation coefficients p,; include the indirect 
dependence of both x; and x; on the n-2 remaining variables. Unless these correla- 
tions are extremely small, the linear dependence of xj and xj is not obvious. The 
direct linear dependence of these two variables is given by the partial correlation 
coefficient 


pil (l,...n)/ij 


This is obtained from the normal equation matrix B (the inverse of the covariance 
matrix) according to 


pig! (1, « » . m)/ij= ~ by/ (dys + By)!” 


Correlation coefficients of varying orders may also be defined, but are beyond the 
scope of this section. 


It is of interest to consider the conditional variance defined as 


G1(t,... mi Vb: 
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in order that it might be distinguished from the marginal variance obtained from the 
covariance matrix M, and used in evaluating the multiple correlation coefficient Rj. 


Consider the bivariate distribution comprised of the variables Xj, Xj- Suppose 
that the error ellipse of a point O is represented by Figure 5-1. oj, oj, are the 
projection of the ellipse on the axes xj, Xj and are the marginal standard deviations, 
whereas the lengths OP, OQ are the conditional standard deviations o (x, x,) and 
o (Xp X;). 


The conditional variance is always less than the marginal variance, unless the 
correlation between xj and x; is zero for all correlation coefficient Rj, defined by 


espe ee ee ype 
oj 


which is a measure of the total linear dependence of the variable xj on the 
remainder. 


5.4 APPLICATIONS AND ANALYSIS 


The preliminary requirement for the performance of the system analysis, is that 
exact data be generated suitable for the combinations listedunder assumption 5. These 
data are readily furnished by the computational sequence listed in Appendix A, 
Construction of Fictitious Model. 


It is to be noted that these calculated data are exact, save for truncation errors. 
As indicated in the descriptive section, the resulting photo-coordinates are not subject 
to displacements caused by aberration, atmospheric refraction, lens distortions and 
the like. This is unimportant, since in practice such displacements would be removed 
from the measured image coordinates, as far as possible, before being subjected to 
the main data reduction scheme. It is however, an extremely simple modification to 
include these displacements should they be required. 


From these computed data one obtains a series of exact photo-coordinates, related 
to the corresponding exact ground coordinates with a transformation involving exact 
flight and camera parameters. The objective of the subsequent computations is to 
determine the effects of various inaccuracies in these coordinates and parameters. 


To illustrate the method of achieving this end, a consideration of the effects of 
errors in the measured photo-coordinates on the determination of ground points 
coordinates, will be made in some detail. 


These are four distinct cases listed under assumption 5. Consequently, there 
will be four sets of parallel computations. To determine the independent effects of 
errors in the measured photo-coordinates on the end results, the influences of camera 
and vehicle path parameters are excluded by assinging to each of them, a zero variance. 








Fig. 5-1 — Error ellipsoid 
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The measured photo-coordinates are assigned variances (and covariances) proportional 
to the specified measuring errors under investigation. Since the system parameters 
are exactly known, together with the transformation relating ground and photo- 
coordinates, the exact ground control data will be assigned zero weights. In this 

way, the resulting variance-covariance matrices of the computed ground points 

reflects only the influence of the errors in the measured image coordinates. 


It is suggested that five solutions be obtained for each caSe corresponding to the 
selected standard deviations in image meaSurement of +1, +2, +4, +8, and +16 
microns. These standard deviations are not quite arbitrarily selected, since they 
correspond to the quoted standard deviations of specific instruments covering the 
range of envisioned mensuration equipment. 


The variance-covariance matrices of the results, is the important output. They 
indicate the effects of the estimated measuring errors or the ground coordinate 
determination. Recalling the previous section, the diagonal elements of these matrices 
represent the marginal variances of the ground coordinates. 


In the event that the computational procedure preserves the normal equation 
matrices, conditional variances, together with the multiple and partial correlation 
coefficients are determinate. However, it is not usual to preserve the normal equations 
in a solution by electronic computer, owing to storage requirements. In this event, 
one can only manipulate with the variance-covariance matrix to furnish correlation 
coefficients, which indicate the effect of all other parameters and variables on that 
one which is being investigated. 


However, remembering that the partial correlation coefficients reflect the effect 
of a Single parameter on the end result, and that the multiple correlation coefficients 
indicate the effect of some selected parameters or variables, one realizes that this 
is exactly what the designed computational procedure accomplishes. Admittedly, if the 
effects of photographic x-coordinates are to be divorced from those of the associated 
y-coordinates, one should assign zero variance to that group of coordinates whose 
influence is to be deleted. This adds two more computational sequences—yielding two 
partial (x, y) correlation coefficients and one multiple (x y) correlation coefficient. 

At present, it is not at all apparent what advantages such a refinement would yield, 
apart from a more complete analysis. This has been mentioned, however, to point 
out that any specified degree of refinement is easily attained. This provides a method 
whereby the effects of measuring errors in the panoramic and frame imagery can be 
separated, leading to valid selections of the appropriate image measuring engines. 
Furthermore, it indicates that retention of the normal equation matrix is unnecessary 
for this portion of the system analysis, implying storage for a larger volume of 
pertinent data. 


At the end of the computational Sequence, one has amassed a sequence of variance- 
covariance matrices, and their normalized equivalents, the correlation matrices. 
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Such numerical data cannot be readily interpreted and evaluated, even by those to whom 
numerical analysis is second nature. 


It would appear that the most striking, and clearly understood, manner of pre- 
senting these individual and collective data, is in graphical form. This has the added 
advantage of being concise and indicative of subtle nuances. The suggested basic 
presentation illustrated by Figure 5-2, is to present a series of graphs whereby the 
abscissa represents the standard error of the photo-coordinates, and the ordinate that 
of the ground coordinates. Correlations, linear, partial, or multiple, are readily 
presented in the form of regression graphs. 


The preceding discussion is readily extended to account for the individual effects 
of each variable or parameter,* so that the basic portion of the system error analysis 
may be considered virtually complete, once the sequence of parameter variation has 
been specified. 


This analysis, however, is restricted to a consideration of those factors explicity 
contained by the transformation formulae. There are certain other factors influencing 
the results, which one has the tendency to ignore owing to their undefined effects. Of 
these, perhaps one of the most significant, is the geometric location of control and 
ground points with respect to camera axes. It would appear that the best manner of 
performing an analysis of this factor is through a judicious selection of control data 
and of points to be determined. The results could be suitably presented in the form 
of a double entry graph—that is with respect to both input variances and angular 
location of the point. 


A similar treatment of the various combinations of ground control data, with 
respect to both quality and location, will also be performed. 


This preliminary error analysis, contained in Table 5-1, is now formulated in 
some detail, according to the following scheme. It is pointed out that the scheme 
Should be applied to each of the four cases listed under assumption 5, and that it is 
easily simplified, or modified to include a more complex model. 


The numerical data obtained from these calculations and their graphical presenta - 
tion, provide a sufficiency of data for a thorough system analysis. Furthermore, these 
data are sufficiently extensive to provide the basis for construction of elegant 


*Henceforth, parameter is used in the sense of a parameter or variable, unless 
otherwise indicated. It can be numerically demonstrated that the results of a con- 
strained adjustment are the same whether the unknowns are considered to be 
parameters or variables. The formal proof of this, has yet to be derived; numerical 
demonstrations indicate that mathematical induction might be the easiest method 
of deviation. 
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Fig. 5-2 — Standard errors oy in position (1), ground X (2), 
Y (3), and Z (4) coordinates as functions of standard error 
in photo coordinates (c,) 








Table 5-1 — Schematic Formulation 


Varied Parameters 


Photo-coordinates 


(xy yy), (x’, y’) 
(xp ya), (x”, y”) 


Inner orientation elements 


(Xpy Yp» £), &p, yp, £), (xp, yp, £”) 


Angular orientation elements 


Sag ase 
Wi, Wi, Wi 


Angular orientation elements 
Ci Vp VF 

Angular orientation elements 
Ki, Kis KF 


Vehicle path elements 


2, w, I, e, T, t 


Pan camera timing 
te, tj 

Ground control quality 
Xi, Yi, 2} 


Ground control quality 
(combinations) 


Xiy Yi, Zi 





Remarks 


Vary individually, in combinations according 
to desired results. Use sequential standard 
deviations of 1, 2, 4, 8, and 16 microns. 


Vary individually, in combinations according 
to desired results. Use sequential standard 
deviations of 1, 2, 4, 8, 16, and 32 microns. 


Vary individually, in combinations according 
to desired results. Use standard deviations 
of 1”, 5", 10%,.30%, 1,2’. 


See 3. 


See 3. 


Vary individually, and in combinations. Use 
a range of variances to cover all expected 
cases. 


Vary individually and in combination, using 
a range of variances to cover expected errors. 


Use 1 all grade I control 

Use 2 all grade II control 
Use 3 all grade II control 
Use 4 all grade IV control 


Judiciously select a reasonable sample from 
all combinations. 
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Table 5-1 — Schematic Formulation (Cont.) 


Varied Parameters 
10. Ground control location 
Xj, Yy Zi 


11. Ground control location and 
quality combinations 


12. Photogrammetric parameters 
and ground control 


13. All system parameters 
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Remarks 


Judiciously select a reasonable sample of 
possible control configurations. Use 
homogeneous data, in each of four sequences. 


Use samples 9 and 10 to construct a 
resultant sample for analysis. 


Assign zero variances to vehicle path 
elements. Use reasonable expected values 
for errors in other parameters, and resonable 
control distribution. 


Assign reasonable expected variances to 
all system parameters. 








prediction and control nomograms. Although the basic system error analysis has 
now been adequately defined, it is thought that the ever recurrent theme of relative 
v absolute accuracy should be discussed. This is the topic of the succeeding section. 


5.5 RELATIVE AND ABSOLUTE ACCURACY CONSIDERATIONS 


It is frequently noted that the use of rigorous error propagation techniques using 
closely estimated variances leads to surprisingly large system errors. As an example, 
if the well known error equations derived by Brandenbeger are applied to a convergent 
strip camera system, which has an extremely high resolution, it appears as if precise 
mapping from these records is not possible. This is a consequence of absolute position 
and orientation errors of the complete system being quadratically combined, so that 
the variance that must be assigned to a small just detectable length may be ten or 
even one hundred times the length itself.Although it is possible to detect and measure 
the images of small objects with great precision, it is difficult to maintain any degree 
of accuracy. This difference, between relative and absolute accuracies, epitomizes 
the different usages of imagery by interpreters and cartographers. In order that one 
might determine this relative accuracy it will be necessary to consider every 
parameter to be exact, (in a two photo-mode’, save five elements which are necessary 
to determine the relative orientation between two pictures. Consequently, ina n 
photosystem it will be necessary to maintain exact values for all but 3n - 1 parameters. 
For convenience, it is recommended that the angular orientation elements, wj, 9}, 

Ky be those that are varied—being virtually independent of any other parameters. 
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Appendix A 
CONSTRUCTION OF FICTITIOUS MODEL 


1.0 INTRODUCTION 


This aspect of the problem is to compute exact data by means of which the 
methods of data reduction and analysis might be tested, and computer programs 
optimized. The construction of these data falls naturally into four sections. 


2.0 DETERMINATION OF SATELLITE POSITION 
AT ANY TIME t; 


For convenience times will be assumed to be given in Greenwich Apparent 
Sidereal Time. Although the time measurement may be in some other system, 
conversions are simple, but necessary. 


Given Data 


Orbital Parameters: 2, w, Ll e, n, 7 
Geophysical Constants: pu = GM = 3.98603 x 107° cm’/sec? (Kaula 1961) 
Observed Time: t; (in GAST) 


Computation 


1. Determine semi-major axis, a, of orbit 


1/3 
a= (u n-’) (111) 


2. Determine mean anomaly, M;, at time t,: 


a. Compute 

yy 27 aT (112) 
b. Compute 

M, = nti + xj (113) 


Po el 
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3. Determine eccentric anomaly, Ej, at time t; 
Ej = M, + e sin Ej (114) 
This is performed iteratively according to 
Mj = Ej (115) 
which is substituted into 


oo ° . ° 
Ej = M; + e sin Ej 


4. Determine the Geocentric Terrestrial Coordinates of the Satellite position 


Xj 


as 
i) 
a 


, according to 


cost; Sin tj0 ][cos® -cosI sinQ [cos w -sin *) -e + cos(E;) 
= 1/2 
r, =a|-sintjcostj0|/sinQ cosI cosQ||sin w cos w| (1-e”) sin E, | (116) 


0 oO I1jL0 sinI 


5. Determine the heading, geocentric latitude, and longitude of the nadir 


-1/2 


COS Ag = Xj (K} + Yi) (117) 
1/2 
sin Ay = Y; (Xf + Yj) (118) 
gt gene Sen Sie age 
cos $y = (x; + Yj) 0G aa ae. Z;) (119) 
2 v2, ot 
sin, = Z, (Xi + Yj + 23) (120) 
cos Ay= sinI cos (Ay - 2) (121) 
sin Ay = coSI/cos $o (122) 
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3.0 DETERMINATION OF CAMERA DATA 


Given Data: Camera half angle a/2 
Ellipsoidal semi-major axis 
Ellipsoidal eccentricity e 
Flying height H 
Camera roll, pitch, wj, 9j 


Preliminary Computation - frame cameras 
1. Determine flying height H 
H=Irj!-R (123) 
(R = radius of earth) 


2. Calculate the azimuth to each of the four corners of the imaged terrain from 


the nadir point Ao, 
Ao, = Ao+ 511/4 = a, Ao, 
Ao, = Ag+ 70/4 = a " 
Ao, =Aj+ I1/4= ag ‘ (124) 
0 


Ay, = Ay + 3 11/4 = or 
Ag, 
3. Compute geodetic latitude of nadir, %), from geocentric $, 
2-1/2 
sin %) = tan 9[ tan? 6) + (1 - e) | (125) 


4. Compute radii of curvature at nadir 


Ny =a (1 - e’ sin? W))-'” (126) 


Ry) =a(1 - e”)(1 - e? sin? ¥)73” (127) 
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Compute mean radii of curvature 


N 
= = 1 
Raja Has Ry sin‘ a, + Ng cos’ a, 28) 
and 
Ra, = Ra = N (129) 


Ry Sin” a + No COs" a 
Compute the angle 6 7 subtended by the arc from the nadir to four corners 


according to 


R, 


8 a = sin! aaa sae -a/2 (130) 
i 


where a/2 is the diagonal half angle of the photography. 
Compute the arc length Sai according to 
Sa; = Ra, Do.» 6a; in radians) (131) 


Compute geodetic coordinates of the four corners according to Rainsford’s 
modification of Clarke’s formula: 


Sq, COS a; 
) % =v) + ——— (132) 
ap mo = Yo + “Ry sin 1” 
b) ¥ = Sa me (133) 
c) compute Ry: 

Ry =a(1 - e”)(1 - e? sin? w)"*”? (134) 
d) compute 
Soy Sin Qj COS a} 
Pi= 2R, Ny sin 1” (135) 








e) compute 


qj = Pj tan a; tan ¥ (136) 
f) compute Sq, (cos a; - 2p;/3) 
a oo 137 
Y= Vo+ Ry sin 1” qj (137) 
g) compute 
_ ee ere oe 
Ny, =a(1- e’ sin’ v,) (138) 
h) compute 
Sea Sai sin (a, - 1p,/3) | 5) 
i= “0+ | cos (G+ 4,/3)Ny sin 1” 
i 


These values Wj, rj define the limits within which the control data for each frame 
camera are to be determined. 


This sequence of computation is performed for both frame cameras. 
For each frame camera, designate the corner points of coverage aS $j, Aj, 


D3, AZ, 3, AZ, Bi, AGG SY, AY, OF, AZ, G3, AY, GY, AY; as in the following diagram: 


oi AY of $224 O7 
Oi 4 Of 3 A3 Og Ag 





Pn 
= NTROL SYSTEM of 





$3,Ag = D3, As 


BY, AG = By, Ay 


Then the area covered by the two photographs is bounded by the points $,, A; 
$2, ra3 $y, A35 %,, rN. 

Let this area be evenly covered with geodetic control points, comprising I evenly 
spaced rows and J evenly spaced columns. 

Then any control point 43;, AGj has coordinates of 


ij = &1 - [(2; - 1G, - &2)/25] -[ 2, - 1) - &,)/2T] 


+ [@; - 1)(2, - 1M, - 2+ 3 - &,)/415] (140) 


ag = 4 ~ [2 - Oy - A2)/25] - [2 - YQ, - 4/21] 


+ [(2; - 1)Q; ~ 1)Qy - Ag+ Ag - AQ)/40T] (141) 


4.0 SPACING OF CONTROL 


Although it is not essential for a super abundance of fictitious data to be generated, 
it is desirable from the point of having superfluous data for statistical purposes. 


It has been stated that each of the panoramic photographs (models) may be divided 
into three abutting segments, each of which will contain a minimum of 5 control points. 
It is considered desirable that there be three columns of control data, one central 
column and two lateral, according to the following diagram. Furthermore, there will 
be seven panoramic models per frame model, which are assumed to be butt joined, 
dictating that 21 columns of control per overlap be established. 


Let it be assumed that the panoramic photography has been rectified to fit the 
frame photographs. If the format size of the frame is a x a mm’, with p-percent 
overlap, then the width of the pan model is 


fig 
100 


‘amm 


~]|— 
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The spacing between each column of control data is 
vB ea mm 


2100 


and the total number of control points per row is 


J 


100 p- 100 
ca rape} * 2 


Assuming p = 60 percent, J = 44 


The interval between each control point is approximately 6 to 7 millimeters, for 
a 9 by 8 inch format. 


Arbitrarily, select the rows of control data 10 millimeters apart, to yield I = 23 
for a 9 by 9 inch format, i-e., 


This leads to 1012 points for the covered area. 


5.0 COMPUTATION OF PHOTO COORDINATES 


5.1 Frame Cameras 


1. Compute $jj, 4;; for each control point within the limits according to 
Equations (141) and (142) 


2. Generate Ix J random integers between 0 and 100 meters. Associate these 
with the $;,, ij in Sequence, to represent the point elevation, Hij, above the 
reference ellipsoid 


3. For the frame exposure times t,, tz, in Greenwich apparent sidereal time, 
compute the geocentric sidereal coordinates, according to 


xX cos tj-sin tj 0 (N+ H) cos gycosaA 
Y} =/sinty cos tj0 (N+ H) cos g Sind (142) 
Z|;; | 0 0 1] [1+ e)N+H sing 


ij ij 


4. Compute the geocentric sidereal positions of each camera station according to 
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Xo cos 2 - cos IT sin cos w -sin w 
sin w cosw 


1 0 re E- e 
0 (1- =| sin E | (143) 


Yo =alsing cosIcos2 b 


ZoJ ts 0 sin I 


as outlined in paragraph 2, this section 


5. For each exposure station, compute the converted photo-coordinates of each 
point, according to 


my, (Xi - Xp) + my. (Yj - Yo) + mys (Z; - Zo) 
xX: =x, -f 11 “Ai 0 12 \Ti 0 13 \Sj 0 (144) 
1 : a Xj - Xp + M939 i7 Yo + M¢33 Zi - Zo 
and 
Mo; (Xj - Xo) + May (Yj - Yo) + meas (Zj - Zo) 
¥,=y- See (145) 
M3 (X = Xp) + M39 (Y; = Yo) + M33 (Z, - Zo) 
in which the values m,, are given in Appendix B, detailed formulation, 
Equations (159) through (167) 
6. Reject all points which fall outside the desired format limits 
7. It is noted that these are corrected photo-coordinates. If desired, apply 
displacements due to aberration, lens distortion, atmosphere refraction, etc. 
5.2 Pan Cameras 


Consider the equivalent frame camera coordinates given by (145), (144) in the 
previous Section, rewritten as 


U; | 
x, = -f WwW, (146) 
and 
= -f ee 147 
ir WW. (147) 
1 
noting that 
Xo = Yo = 0 


eeprecnnerinacine 
er rer Howe ig “ose = ee 
vets < oe foe Fees 
~~TPOL SYSTEM ENTLY cae wang ‘ Sten 2 








Let the equivalent frame be exposed at the mid-point of the scan, at toj. Any 
control point ij; Ay is scanned at some other time t;. By selecting an approximate 
value of tj = tj, approximate values of 4x; = x 


y= 3; 
U; =U; 
Vi = V; 
W, = W; 


may be determined. Consequently one may write 


ais ea a(x;)/at 
o | Yi 
as | aly, )/at 


Ui 


i i - a(U;)/at - etc. 


Since y, is the most rapidly changing function with respect to time, use 
i 


fall 
ti = ti - aye)/at (148) 
as the control equation. This is rewritten as 
=¢° . ° ° ° ° -1 i 
t= tp - (V7/Wi(0(VS/w?)/at) (149) 


as M9 


a as 
pial 


il? Seoaet 








in which 





a a 2m 
— (, - Xp) + (¥, - Yo) ore 8 (Z; - Zo) 


aXi oY; 9Zj 
+ Mo, at + Mg at 4 M93 et. 


ax ay aZ 
o/xx7°) ° 2 pies eat pote £0 
a(V;/W;) = 1/W; Ma >; 7 M2 op 7 Ms Se (150) 





yf | OS Om; - Xo) + SE cw - Ya) + SR (2, - 2) 


aZ; 















aXi aX, eYi  aY i aZ 
pees ee ae aan eet os weg 
+Ms1\ st ot |*™82) at ~ ot | +m at at 
In this equation 
aXi ee . 
ee [(cos w cos Q - sin Q sin w cos I) sin (E - e) + 
(151) 
_ e2)t/2 : Sete 
(1 - e?)!/? (cos 2 sin w + sin 2 cos w cos I) cos E] ieee 
sp [(cos w cos 2 + cos 2 sin w cos I) sin (E - e) + 
(152) 
1 - e2)!/2 (5; ‘ 2 ee. eee 
(1 - e”)'/? (sin Q sin w - cos 2 cos w cos I) cos E] icc er 
“i - a [(sin I sin w sin {E - e) - 
~ elit? n (153) 
(1 - e)!/* sin I cos w cos E] focare 
oXj , 
oa -X, sintj - Y, cos t; (154) 
aY; 
2S ae L 
eer X; cos t; - Y; Sin t (155) 
aZ; 
Fi = 0 (156) 
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The various values of am,)/at are obtained as 








om : am 80; om OK. 
am] _ OMKL Oy KF oi (157) 
at aw at a0; at aK; at 


Since the orientation angles are expressed in the form of 


?; = Por olt; - to) (158) 


a ea Ky + K (t; - to) 
the terms aw,/at, ay,/at, ax ,/at are evaluated as in w, y, K, respectively. 
The expressions dmx /2 Qj, % = Wj, Yj, Or Kj are given by Equations (207) through 
(233), Appendix B, which require the appropriate values of Wj, Pj, Kj as given by 


Equation (158), above. 
The specific computational procedure is: 


1, Assume tj (¢ toi) 
2. Determine X;, Yj, Z; according to (142), for t; 
3. Compute Xo, Yo, Zp at time t; according to Equation (143) 


4. Evaluate m,), according to Equations (98) through (108), Section 4, using 
a: = Ajo5+ @ (ty - toj) 


J 
5. Compute V?, Wj and ¥; , according to Equation (147) 
6. Compute aX;/at,.... aZ/at, according to Equations (151) through (156) 
7. Compute am,,/a a; , according to Equations (207) through (233), Appendix B 
8. Compute the six equations am,,/at according to Equation (156) 
9. Compute a(V;/W;)/at according to Equation (150) 


10. Evaluate Equation (149): 


tr = tp - [ acvi/w))/at ]-* (vi/w)) 


1 1 
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11. Put t°° = t? 
i i 
12. Repeat steps 2... 11 until the value of 
° ° ° ° -1 
(Vs/W5) (2 (Ws/W5)/at) 


is negligibly small. 


This is done for each of the points on each of the equivalent frame pictures. These 
points are then transformed into the appropriate coordinates x.’, Yj , according 
to Eqs. (19) and (20) in Section 2. 


1 








Appendix B 
DETAILED FORMULATION 


The terms of the orientation matrix M = MjOj, are suppressing the subscript j on 
all except the rotation Wj. 


™Mi1 


-cos Ay {[cos k cos y cos Ag+ (sink cos wj + cosk sin gy sin w;) sin Aol 
sin ) - (sink sin w; - cos x sin y cos 4) cos 6} 
+ [cos k cos ¢ sin Ay + (sin x cos wy + cos x sin g sin w) cos Ao] sin Ap 


(159) 


-COS Ay {[-sin K cOS g cos Ag+ (cos x cos wy - sin k sin g sin wj) sin Aol 

sin - (cos x sin wj + sink sin g cos ;) cos &} 

+ sin ) [sin « cos g sin Ay + (cos x cos wj - sin « sin 9 sin w;)cos Ao] 
(160) 


“COS Ay {[sin 9 cos Ao - cos ¢ Sin j sin Ao] sin & - cos ¢ cos wj cos d} 


+sin Ay [-sin ¢ sin Ay - cos ¢ sin wj cos Ao] (161) 


-Sin Ag {[cos K COS ~ COS Ag + (sin k cos wj + COs kK Sin g sin w;) sin Ao] 
sin , - (sin « sin wj - cos « sin g cos j) cos Sy} 
-cos Ay [-cos k cos w sin Ay + (sink cos wj + cosx sin y sin w;) cos Ag] 


(162) 


-sin A, {[-sin k cos @ cos Ay + (cos Kk cos wj - sink sin ; sin v) sin Ag] 


sin $, - (cos x sin wj + sink sin y cos w;) cos bo} 


-cos Ay [sin k cos @ Sin Ay + (cos k cos w) - sink sin @ sin w,) cos Ao] 


(163) 








M3, = -Sin Ap [(sin gy cos Ay - cos ¢ sin wj sin Ay) sin &) - cos g cos aw) cos $y] 


-COS Ag [- sin gy sin Ay - cos ¢ sin Wj COs Adl ‘ (164) 


My = [cos x cos @ cos Ay + (sin k cos w;+ COs K sin gy sin w;) sin Ay] cos #, 


+ (sin « sin wj - cos k Sin y cos w,) sin (165) 


m3 = [-sin k cos y cos Ay + (cos K cos w; - sin x sin g sin w;) sin Ao] cos , 


+ (cos k sin w; + sin x sin ¢ cos wj) sin $, (166) 
m3; = (sin cos Ag - cos ¢ Sin j sin Ay) cos $y + cos g cos wj sin &, (167) 
In which 
COS Ap = Xp/ (XK? + ¥3)!/2 cos &y = (X4+ YX + v2 4 ZA) 
sin Ay = Y)o/(x? + ¥3)!/2 sin ) = ZX} + ¥3+ 22)” (168) 
cos Ay = sin I cos (Ay - 2) sin Ay = cos I/ cos % 
OF ,/dx; =1 aF,/ax; = 0 
aF,/ay; = 0 aF,/dyj = 1 
OF,/8x, = -1 aF,/ax, = 0 (169) 
aF,/dyy = 0 8F,/dyo = -1 
OF\/ef = U;/Wj OF,/af = V4/W; 
where 
U; = mys (Xj - Xo) + myy (CY; - Yo) + myg (Zj - Zo) (170) 
Vj = Mgy (Ky - Xq) + gq (XY; - Yo) + mg3 (Z; - Zo) (171) 
Wi= M3; (Xj = Xy) + Myo (Y; = Yo) + M33 (Zi = Zo) (172) 


For pan camera it is noted that it is more appropriate to use the expressions 


U; 
F, = tan aa 








and 


, Vi 
F,=(y -y,)+f Wi 


whence one obtains 


aF,/ax; = - sec” aj/f’ 


aF,/ax, = + sec? aj/f’ 


OF ,/ayy = 0 
aF,/af’ = - a; sec? a;/f’ 


Note that although 


C= 


a 
RQ | 


aF,/ax; = - C [1 + ly’ - yg) tan? aj]/f’ 

aF,/ax, = + C [1 +(y’ - yg) tan? a; )/f’ 

aF,/ay; = + sec aj (173) 
OF ,/ayp = - sec Qj 

aF,/at’ =(v/w) - aj [C - ly’ - yg) sin a;]/f’cos? a; 


C is not a variable in the expression F,, and must be considered invariant. Further- 
more, it is important to note that for a pan camera, the factor (f’) assuming the 
subsequent partial derivatives should be made equal to unity. 


Working in sidereal system, i.e., put 


Xj = Xj cos tj - Yj sintj, Yjg = Xj Sin ty + Yj cos ty , Zig = Zj 


ae = = [ms cos tj + my sin i - ay (mas cos tj + ms, sin t)| (174) 
so = = [-mis sin ty + my cos tj a (-mey sin tj + mz cos t)| (175) 
ae ~ : [mis "Ww ms | (176) 
= == [+m sin t; + mz Sin tj = (+m 3; cos tj + mgz sin t)| (177) 








oP eat [-mas sin tj + my, cos t; - 3 (omg, sin t, + ms; cos t))| (178) 

















9Yis Ww 
oF, _ tf V (179) 
dZig W (me W may 
Define 
xX Ay + By is (E) - i 
¥| =|Ay + By| [sin E = Ry (180) 
Z Az, +B, 
Ay = a(cos 2 cos w - sin Q sin w cos 1) 
Ay = a(sin 2 cos w + cos Q sin w cos I) 
A, = alsin I sin w) 
& is (181) 
By = -a(l-e?)"* (cos Q sin w + sin 2 cos w cos I) 
By = -a(1-e2)'” (sin Q sin w - cos 2 cos w cos I) 
Bz = a(1-e*)”” (sin I cos w) 
ax® . : 
Tea {(sin 2 cos w + cos 2 sin w cos D[cos (E) - e] 
- (1-e2)'” (sin @ sin w - cos 2 cos w - cos 2 cos w cos I) sin E} (182) 
ay? z : 
sate {(-cos 2 cos w + sin Q sin w cos I) [cos (E) - e] 
+ (1-e?)!” (cos Q sin w + sin Q cos w cos J) sin E} (183) 
az? 
sq 8 (184) 
ax? : - 
35 772 {(cos Q sin w + sin 2 cos w cos I) [cos (E) - e] 
+ (1-e2)'”? (cos 2 cos w - sin Q sin w cos I) sin E} (185) 
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where 








i ee {(sin 2 sin w - cos 2 cos w cos I) [cos (E) - e] 


aw 
+ (1-e2)” (sin 2 cos w cos 2 sin w cos I) sin E} (186) 
: 1/2 
- = +a {sin I cos w [cos (E) - e] - sin I sin w sin E (1-e?) pa (187) 





= = a{sin Q sin w sin I [cos (E) - e] + (1-¢2)” sin 2 cos w sin I sin E} (188) 








~ =a {-cos Q sin w sin I [cos (E) - e] - (1-e”)”? cos 2 cos wsinI sinE} (189) 
aZ° . 9\1/2 ; 
or a{cos I sin w [cos (E) - e] + (1-e?)" cos I cos w sin E} (190) 


at =a{(cos 2 cos w - sin 2 sin w cos I) [Ae] 


- [Be] E (cos & sin w + sin 2 cos w cos D} (191) 


a = a{(cos 2 cos w + sin 2 sin w cos I) [Ae] 


- [Be] (sin sin w - cos 2 cos w cos I)} (192) 
aZo . : WN [ ] ; ( 
36 72 isinI sin w e] + [Be] sin I cos w} 193) 
[A ] Son sin? E 

€1* "lt F1-e cos E 


_ sin E (cos E - e) 
[Be] = (i-e cos E)(1-e?)!? 


It is noted that the problem will be concerned with observed times t;. Since M isa 
function of t;, the time of epoch 7 will be usedasan independent variable. 





cf — 

















ax? 
= = , -a[(cos » cos Q - sin Q sin w cos I) sin (E) 
ti -7 
_a2)1/2 : : i 
+ (1-e7)! (cos Q sin w + sin 2 cos w cos I) cos E] (cos ED 
-2 Rx pis 1 °3/8a (194) 
oy =, -al(sin Q cos w + cos Q sin w cos D) sin (E) 
(1-e2)1/2 (53 : Z ie 
+ (1-e2)! (sin 2 sin w - cos 2 cos & cos I) cos E] pt 
-2 Ry 118 7°°/8/3a (195) 
eae -a [sin I sin w sin (E) - (1-e*)!”” sin I cos w cos E] al 
an (i-e cos E) 
-2 Ry ul3 75/3/30 (196) 
ax° ex® - 
~ OX, _=7_ (197) 




















a7 tg T age) 
aZ 22? =n | 
Ta eee (199) 
Fi my 3 , 
il +f omy a : am am 2 
a0; Wi {0% Sa Ta Boy] > Yo aay ay 
Q 
+(Z,-Z off m8 Ty al (200) 
“7 j 
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RET 











Fi +f E am | E am 
2 - Ty; —it] 4 (y, - Y,) | —# - r,, — 
ag 3 ox, - x0) ag Hag | TMi TY Lag “TH ap 
3 
+ (Zj - Zo) ES Ty omsl| (201) 
OF, +f dm am am am 
Fi tf Vy. - op. Mal iy, —y,) {2m p, Ome 
ow fox x,) | Dy sae] + v0) mi Ts | 
() 
+ (Zs - Zo) sz - i; “pa]| (202) 


OF; re) ) 
Foi if 2M, a [2m Ty oma] + (yy -¥) E: Ty oa 
ej J 


aw; Wj 9Wj ow; 
() 
+ (2 ~ 2) [2B - Tas al (203) 
J 


oF.j _ tf omy dm dm 
—2U = toe, - 20 [3 Pe re 3 2miat) «(xy ~ ¥o) [Be - ry PH] 


ap Wi ag ap 
8m dm 
‘ ci: ig ree 
+(Z, - Zo) E Ty ae } (204) 
oF,j +f am am om am 
pele Seyi ar 3 ads goto ie am .r,, ia 
= i lex, Xp) | Tae say + (Yj; vo) [° lai | 
am am 
+(Z, - Zo) [zee - Ty ‘za! (205) 


or in general 


SF eg, - x) [BL ry 2B] + xy = ve) PERRY ry 2 


ae Wy FY FY ae 
+(Z, = Zo) Ee - Ty al (206) 
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where 
r., = @ki (Xj - Xp) + mp (Yj - Yo) + tg (Zi - Zo) 
ki © mg, (Xj - Xo) + mgy (Yq - Yo) + mg3 (Zy - Zo 
k= 1,2 


am ; ; 
aa cos Ap [(sink cos wj + cos x sin @ sin w;) cos %p 
J 

- (cos x sin y cos w - sin x sin w) sin Ay sin %9] 


+ sin A, cos Ay (cosk sing cosw-sinkxsinw)  . (207) 


dm ; : ‘ ; 
oo sin Ay [(sin k cos w + cos x Sin g sin w) cos $y 
J 
- (cos « sin @ cos w - sink sin w) sin Ap sin ®o] 


- COS Ay COS Ay (cos k Sin g cos w - sin k sin w) (208) 


smd = (cos x sin g sin w + sin x cos w) sin $y 
j 


+ (cos x sin yg cos w - sin x sin w) sin Ay cos (209) 


om ‘ ‘ ; : ; 

ae cos Ay (cos x sin w + sin x sin y cos w) sin Ay sin %, 
J 

+ (cosx cos w - sin x sin g sin w) cos ®o] 


- Sin ») cos Ay (cosk sin # + sin K sin @ cos w) , (210) 


am : : i i i i 
smut = sin \y [(cos « sin w + sin x sin y cos w) sin Ay sin ®y 
j 


+ (cos x cos w - sinx sin @ sin w) cos #5] 


+ cos \(cosx sin w + sink sin @ cos w) (211) 


dm : ; : 
ea Se eOs x COS W - Sin x Sin @ sin w) sin 
J 


- (cos k sin w + sink sin g cos w) sin Ay cos %) (212) 
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_—10P SCE 
ae 








a mar : 
=i 1 = cos Ay [cos ¢ cos w sin Ay sin %y - cos y Sin w cos 99] 
j 
- sin A, cos g cos Ay 
3 A 
at? = sin Ay [cos ~ cos w sin Ag sin $y - cos g sin w cos ®o] 


+ COS Ay COS ~ COS w cos Ag 


am : : : 
—i = - cos g cos w sin Ay cos ®p - cos g sin w sin $y 


am : : , , 
<Tit = cos Ay {[cos x sin y cos Ay - cos K cos g Sin w sin Ag] sin &p 


-cos kK Cos Y cos w cos $,} 


+ sin Ay [cos x sin y sin Ay + cos x cos ¢ sin w cos Ag] 

a 

i = sin \y {[cos xsin gy cos Ag - cos x cos ¢ Sin w sin Ag] sin %y 
- COS K COS Y cos w cos &p} 


- coS Ay [cosk sin y sin Ay + coSk cos ¢ Sin w cos Ag] 


) ‘ Z ; 
—13 = cos $, (cos «cos y sin w sin Ay - cos xk sin Y cos Ap) 


- cOS K COS Y cos w Sin 


= cos A, {[sink cos g cos w cos %) 





-[sin x sin.g cos Ay + sink cos ¢ sin w sin Aj] sin ®p} 


- sin [sin x sin @ sin Ay + sink cos @ sin w cos Ao] 


(213) 


(214) 


(215) 


(216) 


(217) 


(218) 


(219) 








——22 = sin A, [(sin « sin w cos gy sin Ay - sin k sin gy cos Ay) sin $y 


+ sin x COS Y cos w cos %] 


+ cos A, (sin x sin g sin Ay + sin x cos ¢ sin w cos Ay) (220) 


am : : ‘ ; ; 
—4i = (sinx sin g cos Ay - Sink cos ¢g sin w sin Ag) cos ®p 


+ sin K COS Y COS w Sin %, (221) 


a 
a = sin Xy (sin g sin w cos Ay - cos @ sin Ag) 

- COS Ay [(cos g cos A, + sin g sin w sin A>) sin $5 

+ sin g cos w& cos 8p] (222) 
a 
me COS Ay (cos ¢y sin Ay - Sin g sin w cos Ap) 

- sin Xp [(cos g cos Ay + sin y sin w sin Ao) sin 8) 

+ sin g cos w cos ®)] (223) 
a 
a (cos g cos Aj + sin g Sin w sin Ay) cos ®y 


- sin Y cos w sin &, (224) 


am : ; ‘ 
es cos Ap {[sin k cos y cos Ay + (sink sin @ sin w - cos K cos w) sin Ao| 


sin @,+ (cos x sinw + sink sin g cos w) cos @o} 
+ sin Ay [sin « cos y sin Ay + (cosk cos w - sink sin g sin w) cos Ag] 


(225) 








am F ; : : ; : 
se = sin Ap {[sink cos @ cos Ap + (sin x sin ¢ sin w - cos k cos w) sin Ag] 


sin 6, + (cos k sinw + sink sin w cos w) cos &p} 
- cos \» [Sin k cos y sin Ay + (cos x cos w - Sin k Sin © sin w) cos Ao] 


(226) 
om : ; ’ ; : 
—13 =[(cos kK cos w - sink sin @ sin w) sin Ay - Sin K CoS @ cos Ag] cos % 


OK 


+ [cos x sin w + sink sin g cos w]| sin %) (227) 


3 P { : ; 
ute cos A» {[(cos k cos y cos Ay) + (sink cos w + cos x sin g sin w) sin Ao] 


sin + [cos « sin g cos w - sink sin w] cos ®o} 


+ sin Ay [cos kK cos @ sin Ay - (sink cos w + cos xk sin g sin w) cos Ag] 








(228) 
F) 
ae = sin Ay {[cos k cos gy cos Ay + (sink cos w + cos K sin g sin w) sin Ag] 
sin 4, + (cos Kk sin g cos w - sink sin w) cos &} 
- cos Ay [cos K cos g sin Ay - (sink cos w+ cos« sin gy sin w) cos Ag] 
(229) 
F) 
ae = cos ®, [-cos K cos y cos Ag - (sink cos w + cos sing sin w) sin Ao] 
+ (cos k sin 9 cos w - sink sin w) sin 4, (230) 
om3; _ 
se = 0 (231) 
am 
sk = 
ee (232) 
am 
——-s2 = 
aK 0 (233) 
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{- [cos k cos y cos Ay + (sin k cos wj + cos kK sin g sin wj) sin Ap]sin % 


+ (sin « sin w; - cos x Sin g cos w;) cos  ,} a(cos Ay)/aN 

+ {-cos K cos g sin Ay + (sin k cos w; + cos k sin g sin wj) COS Ao} * 
8 (sindj/an 

+ {-cos Ay [cos k cos y cos Ay + (sin k cos w;) + cos x sin y sin w;) 
sin Ao|t x 8 (sin &,)/aN 

+ [cos Ap (sin « sin w} - cos k Sin y cos w;)] a (cos &y)/an 

+ {-cos Ay [cos k cos ¢ sin &) + sin Ay 

(sin k cos j + cosx sin g sin wh x 8 (cos Ay)/a2 

+ {-cos Ap [(sin k cos wj + cos x sing sin w;)] 


sin ® - sin 4) cos g cos x} x 2 (sin Ay) /aQ (234) 


{cos kK cos ¢ sin Ay - (sin kK cos wj + cosk sin ¢ sin w,) cos Ao} x 


8 (cos Ay) /ak 

+ {-[cos k cos y cos Ay + (sin kK cos w; + cosx sin g sin w}) sin Ag] 

sin $, +(sin k sin w; - cos kK sin @ cos w;) cos $4} a (sin Ay)/a 

+{-sin Ay [cos « cos y cos Ay + (sin x cos w; + cos sin ¢ sin w,) 

sin Ao]} x 2 (sin &)/a2 + {sin dy (sin « sin wj - cos x sin gy cos w;)} 

8 (cos &))/a2 | 

+{-+sin Ay COS K COS @ Sin Py - COS Ay (sin k cos wj + cosk sino sin w;)] 
x 8 (cos Ay)/a2+ {-sin A» sin & (sin « cos wj + cos x sing sin w)) 


+ cos Ay cos Kk cos gy} 2 (sin A,)/a& (235) 


[cos K cos ~ cos Ay + (sin k cos wj + cos K Sin ¢ sin w) sin Aol 


8 (cos B))/a 
+ [sink sin w; - cos« sin g cos ,] 4 (sin $9)/a2 + [sink cos w; 
+ cos k sin g sin w;| cos $9 * 2 (sin Ay)/22+ [cos cos y cos 4] 


8 (cos Ay)/a2 (236) 


ss, 7 an Sh Bee at 
Sere Sh eae 
a er ae ae ery 





“Tut = {(cos « sin wj+ sink sin g cos w;) cos $,+ [sink cos y cos Ay 


- (cos k cos w; - sin« sin y sin w;) sin Ag] sin ®y} 3 (cos Ay)/aN 
+ [sin x cos g sin Ay + (cosx cos wj - sink sin g sin wj) cos Ag] 
8 (sin Ag) /a2 

+ {cos Xo [sink cos ~ cOS Ay - (cos k cos wj- sin x sin gy sin wj) 
sin Ao]} 2 (sin &))/a2 + {cos ry [sin x cos ¢y sin %p] 

+ sin Ay (cos k cos w; - sink sin gy sin w,)t 8 (cos Ay)/a2 

+ [sin Ay sink cos y - cos Ay (cos cos w; - sink sin g sin w;) 
sin $,] 2 (sin A,)/a2 


+ cos y (cosk sin wj + sink sin g cos wj) 8 (cos ))/a2 (237) 


—2 = {(cos k sin wj + sink sin g cos w;) cos ®, + [sink cos g cos Ay 


- (cos x cos w; - sin « sin w; sin g) sin Ao] sin $9} 2 (sin A9)/an 


- [sin k cos 9 sin Ay + (cos x cos w;- sink sin gsin wj) cos Ag ]2(cos A9//aN 


+ [sin Ay (cos x sin w;+ sin x sin g cos w,)] 8 (cos )/a2 


+ sin A) [sin k cos g cos Ay - (cos kK cos w - sink sin w sin 9’ sin Ag] 
8 (sin &))/a2 + [sin Ay sin x cos —y sin %p - cos A 


(cos k cos w, - sin « sin g sin w,;)] 8 (cos Ay)/an 


- [sin X% (cos x cos w; - sin k sin w; sin vy) sin 


+ cos Ay sin Kk cos g] 9 (sin A,)/an (238) 


dm P , : ; : 
aryl co kK COS w; - sink sin @ sin wj) sin Ay - sin x cos v cos Ay] 


3 (cos $,)/22 + [cos k sin w; + sink sin g cos wi] 8 (sin 6))/aQ 


- [cos ) sink cos g] 2 (cos Ay)/aQ + [(cos « cos w, 


- sin « sin g sin w;) cos 6] a (sin Ay\/an (239) 
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am 
1= [cos y cos w; cos %, - (sin g cos Ay - cos g sin w; sin Ag) sin ®p 
3 J 


9mM32 


J 
2 (cos Ay'/a - [sin g sin Ay + cos y sin w; cos Ag] @ (sin Ay)/aN 
+ [cos X» cos g cos wj] 2 (cos $9)/ak - {cos Ay [sin ~ cos Ay 
- cos ¢ sin wj sin Ao]} 2 (sin @,)/a2 - [sin Xy cos ¢ sin OF 
+ sin @ Cos Ag sin yo] 3 (cos Ay)/8N + [cos Ay cos gy sin wz sin bp 


- sin Ay sin ¢g] 8 (sin Ay)/aQ (240) 


ae [sin g sin Ay + cos g Sin w; cos Ag] 2 (cos Ay) /ak 


J 


+ {cos ~ cos w; COs ®, - [sin g cos Ay - cos ¢ sin Ww; sin Ay)] sin ®} 


9 (sin Ay) /9Q + [sin Ay cos yY cos w;] 8 (cos ®p)/an 
- [sin Ap) (sin g cos Ay - cos @ sin w; sin A,)] 2 (sin ®))/a2 
+ [cos Ay cos g sin w; - Sin Xp sin y sin ®)] 8 (cos Ay)/an 


, sin ®) + cos Ay sin g)] 2 (sin Ay /an (241) 


+ [sin Ay cos ¢ sin w; 


rs) 
= = [(sin g cos Ay - cos g sin w; sin Ay)] 2 (cos &,)/a2 
+ [cos ~ cos 5] 8 (sin 6,'/aQ + [sin v cos ®y] 3 (cos Ag) /a2 
- [cos g sin w; cos #9] 2 (sin Ay)/a2 (242) 
a (cos rg) Y ax ay 
@ \cos Ag? __ YY OX OS OY 
32 (x? 4 ys lx an 59 (243) 
d(sin dr) _ _X — (cos rp) (aaa) 
aQ Y aQ 
a (cos ®)) sin ®, tan & ax oY 
= —_ — 
aes en eae [x 3 ae er (245) 
d (sin $,) a (cos #,) 
pda eis | Laan peered 
a cot dy ——=S (246) 
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8 (cos Ay) _ sinI 7 X sin 2+ Y cos 2 ales) (247) 
ary am a sin (A - 2) Se ge a2 
o Asin Ag (sin Ad) _ _ Sind 2 (cos #) (248) 
a2—é—“‘écc#OS OCA 
@ (cos Ay) _ ¥ ax | ( 
<2 Ey a -xX = 249) 
a(sin ry)  _ X_ 2 (cos Ay) (250) 
aw Y aw 
8 (cos $)) _ 1 ay az 
eat = dr {sin 25 tan & [x % oy 2] - 2 eos #22 (251) 
d(siné,) 1 aZ ss sin ® ax 
gksin @y 5 2 ce) Le BBO pla 
= = {cos Fic =; aut (252) 
a (cos Ao) _ sin 1(X sin @ + ¥ cos @) , 2 (cos Ay) (253) 
aw 7 Y 
(sin Ay) __ sin Ay 2 (cos #9) (954) 
aw cos ®y aw 
8 (cos Ao) _ Y ax . aY 
at + [y Tie | a 
a (sin As) __ X __ 8 (cos Ag) 
al20t:té<‘CY:~C*é«CY uaa! 
a(cosd) 1). ax ay Z 
ae eam RZ isn B tan Bo [x 81 +¥Y =] - Z cos Do 4 (257) 
a (sin é)) 1 eZ sin & ax 
pe cea = 2 —_ hacen} — — 
al 5 feos oor RR Tae ee) 
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a (cos %o) (959) 


8 (cos Ag) _ cos I cos (Ay - 2) + a * (x sin 9 +¥ cos ) al 








at 
d (sin Ao) 1 : ? 3 (cos #0] ( 
eee Se ee = 260) 
al cos &, | Reelee al 
a(cosr))_SsiéY ax . oY 
2 Moos 0 = an [YS x Z| ey) 
a (sin X) _ _X _ 2 (cos Ag) (262) 
ae Y de 
a (cos &,) er Oe : ax 24 P aZ 
era Az {Sin Sp tan Sp |X 5 + ¥ ZU] - Z cos o 5 (263) 
a(sind) _1) 2,92 sin [, ax 
oS R cos’ ® 5, R \* de * ae (264) 
@ (cos Ay) _ gat [X sin 2 + ¥ cos Q] 2 (cos Ay) (268) 
de Y de 
3 (sin Ad) 2. sin Ao 7 ) (cos $y) (266) 
de cos By de 
a (cos Ay) Y aX . oY 
ee vtedeteah aa sent y ye ee 
GSE | oy * ay | ae 
(sin X») _ _ X | 2 (cos Ay) (98a) 
an Y 2am 
a (cos) 1 ax oY aZ 
se ; paints eee oe ae 
op 3 si #, tan [x an +¥Y a] Z cos 8 mn (269) 
a(sin®)) 1 aZ sin ® ax ay 
Yo 4+) 2g, 24 _ SMPo |y 9A na 
aR = fos es = xz + ¥ = | (270) 
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8 (cos Ay) _ — [x sin 2 + ¥ cos 2] 2 (cos Ay) (S94) 


an Y: an 
a(sin Ay) __ sin Ag 2 (cos $4) (979) 
a7) cos $y an 
9 (cos Ag) _ Y ax | 2x] 
ar xe 2) 2 [* a7 Xr (273) 
a (sin Ay) __ X a (cos Ag) (974) 
aT Y aT 
a(cos%) 1), ax ay] _ az 
a7 7 Rr Sin So tan & |X" + ¥ 57] - Z cos & 7 (275) 
d(sind,) 1 2g 22 sind [lax . ay 
or RCS re Boe ae (276) 
@ (cos Ag) _ as [x sin 2 + ¥ cos Q] 2 (cos Ay) (277) 
aT pf oT 
a (sin Ay) __ sin Ay 2 (cos $5) (278) 
OT cos $y aT 


Now prepared to form the partial derivations of F,, F, with respect to the orbital 
elements. For any element &, may write 


8Fy +f E _U W] (279) 
ag 
and 


(280) 


ecb ahbbdne 
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which became 


7 E ™M2 ~ ms an E M33 - mis| ozs (281) 
and 

es oe Eas 7 a] * [z M31 ~ M21 = 

= EF M32 = mz] a¥0.,(E M33 ma| 27a | (282) 


generalized to 
OFei Hf iy _ [Be - _ 9M3y - oma «Ty; 2m] 
ae (X; - Xp) aE Nki : | * (Y, Yo) BE ki “9¢ 


3m dm aX 
+ (Z; = Zo) Ee Vii | + (Tig M3; ~ my) BE 


ay az 
+ (Dygq mg - mye) “5p + (Tig mg ~ my) 22a (283) 


where I; has been defined previously 


Now all terms have been defined — coefficient matrices can now be formed. 








Appendix C 
INVERSION OF LARGE MATRICES 


1.0 INTRODUCTION 


One of the stumbling blocks preventing the rapid adoption of analytical photo- 
grammetry is the problem of solving a large system of linear equations, which is 
intimately related to the problem of inverting a large matrix. 


Consider a system of n linear equations relating n unknown quantities. Let 
the system be represented as 


AX=R (284) 


where A signifies the coefficient matrix 
R signifies the vector of constants 
X signifies the vector of unknowns 


Provided that the matrix A is not singular, the determination of the inverse matrix 
A~! and of the individual unknownsx,;, may be readily obtained in theory. However, 
if A is a matrix of high order, a formidable amount of computing is required to 
obtain the solution. In general, one requires not only a specific solution for the 
vector X, but also a theoretical error analysis of the results. For this, it is 
necessary to determine the inverse matrix A-}. 


With the available high speed digital computers the problem is essentially one 
concerning the storage requirements for the quantities involved in the computations, 
the time taken to perform the calculations, and the precision of the results. This note 
will be an attempt to determine the optimum procedure which will furnish the required 
solution for X and for A“. 


The methods by which the inverse of a matrix may be determined are defined as 
direct if the inverse may be determined by a finite number of operations, and as 
iterative if the inverse is determined as the limit of successive approximations. The 
direct solutions do not yield an exact solution, on account of the accumulated round - 
off errors. These errors may be considerable in the case of a large matrix. It will 
be necessary to determine the accuracy of the determination of A~', and to provide 
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a method of correcting the calculated elements. In view of this, only the direct methods 
of evaluating the inverse will be considered, and a method of correcting the inexact 
solution will be given. 


2.0 HOTELLING’S METHOD OF CORRECTING AN 
APPROXIMATE INVERSE MATRIX 


Let the approximate inverse of the matrix A be D’, and define 


F°=I-AD° (285) 
Note that if 

D°= A! 
then 

F°=0 


Now consider the sequence of matrices 

Dj = Dj-1 (1+ Fy-4) (286) 
and 

Fj = (I - AD,) (287) 
where i= 1, 2, 3,...,m 
so that 

Fr = 1 - AD = F2y-1 = Fy-p =... = FO™ (288) 


Provided that the norm of F° is less than or equal to unity, the sequential determina - 
tions of D, converges rapidly. This is based on the convergence theorems: 


If a system of linear equations is written in the form X = AX° + B, the covergence 
of the iterative solution for the vector X starting with an initial vector X° and with 
any value of B is guaranteed if the proper numbers of the coefficient matrix A are 
less than unity in modulus. 


From this theorem it follows that: 


In order for the iterative process to converge, it is sufficient that the norm of the 
matrix A be less than unity. 


C-2 








3.0 MATRIX INVERSION BY THE SINGLE 
DIVISION SCHEME 


This method is based on the theorem which states that if the leading submatrices 
of the square matrix A are non-zero, then A may be factored into the produce of a 
lower triangular matrix by an upper triangular matrix. 


Consider the matrix A and the upper triangular matrix B and the lower 
triangular matrix C, such that A = CB. 


This matrix C may be formed by using the recurrence formula 
j-1 
Cij = ij - 2s (cays), i> i, Cik = ik (289) 
k=1 


and the matrix B by the recurrence formula 





i-1 
aij - be (ci_by;) 
bij = a Pi, ayed (290) 


Since A~! = B=! C7! then A“! C = B7!, in which B™ is an upper triangular matrix with 
1 along the leading diagonal. Similarly, BA7!=C-!, The equation for B~! furnishes 
n(n+1)/2 equations of the form: 


ei 


ee ee ee : oe 
11 jl] 21 i2 


(8... fice ee. a St 
31°33 nl “in 


-1 


“1 “i _ 
Co2 Aig + Caz Aig t+ ++ Cg antl 
C33 aig teste antl (291) 
“1 _ 
Cnn 4jn = 1 
eed Geen 0 


Similarly, the equation for C™! furnishes n(n-1)/2 equations of the form: 


deniers es es 


s we {iA FACE NT-REY HOt 
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“4 “1 “1 “1 
Ag; * Oya 49y Pug Pa Pin ay 
ai+b +b. at =0 


=i 
aj * 23 433 7° °° 7 Pen Anj 
1 “1 2 
oY ale a a a 0 
(292) 


st ofs 
an-1j 7 ce Ani 


j=2...n 
From these two sets of equations the required values of aj} are determined, by 


putting i =n in Equation (291) and solving for ann . .. az}, which are then substituted 
into Equation (292) to obtain a(}_3), for j =n. 


A symmetric matrix S may be factored into two triangular matrices such that the 
lower triangular matrix is the transpose of the upper triangular matrix, i.e., 


s=D!pD (293) 
Then 
st=p7?pr-l (294) 
According to Equation (290) 


i-1 
Sij - Se (dy idy;) 


1 
dy = —— ij (295) 


ai, = Sij - » (dye)? (296) 








Let D = @B, such that 


os ee 
D= Pr 1 a on (297) 
co 1 
The elements of D-! = di} are 
djj = 0 for i> j (298) 
dif = aa a a (299) 
n 
dij = > bi di; for j>k (300) 
k=i+1 


4.0 PARTITIONING METHOD 
Let the square matrix S, (# 0), be partitioned according to the following scheme: 





aN | : B 
S ca ( ) 
= | p+q=n,p2>q 
Cc q! D 
+ | 


and let the inverse S7! be partitioned correspondingly according to 


M N 


Since S.S7! =I the following relations are evident: 


7 Peep Sree Vere 


~TARENT-Ke (HOLE - 
CONTROL SYSTEM ONLY 








Set I 


N = (D-CA-! B) 


= -1 
M =(-NC A“) (301) 
L =(-A-! BN) 
K = (A7! -A7! B M) 

Set I 
-1 

K = (A-BD™ C) 

. -1 
L=(KBD™) (302) 
M = (D"'C K) 


N = D-! -D-! CL 


The solution involves the inversion of one p < p matrix and one q < q matrix; 
various matrix multiplications. 


5.0 BORDERING TECHNIQUE 


This is based on the partitioning of a n x n matrix A, into a (n-1) x (n-1) 
matrix An., whose inverse is known, a row matrix r, a column matrix c, anda 
single element, apy, according to the following scheme. 


‘ 
An-1 , An-1 1 Cp 
An= ; = . (303) 
| 
jon-1.n l 
see Tee ec 
{ 
The correspondingly partitioned inverse is 
1 
APy-1 Sn 
-1 
An a Sere ee | = (304) 
Gn | i/o, 
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Since 


I 0 
AA™!= 
0 il 
then 
6 ees (305) 
an 
@n = 4pn - Pn An-1 Cn (306) 
“rp An 
Pg ea Hak (308) 
an 
so that 
-1 -1 ~1 
re : An-1 €n Tn An-1 . An-1¢n 
-i = On ee: | a, 
An = | (309) 
At! 
Tn “n-1 1 
a a 


The inverse of the matrix is thus obtained by successive borderings, through the 
successive inversions of the sequence: 


Ay > {811 Ag2] > [Ast Biz 4g 
21 Age Af. 822, Sag | aes 


431 gn «= 3 


Each step of this determination is carried out according to Equation (304), 


In the event that the matrix to be inverted is symmetric, the bordering technique 
becomes much simpler. 
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Since the ith row of the matrix is identical with the jth column, the corrective 
terms 


Aqt-1 Cn-1 Mn-1 An-1 
an 


which are applied to Ant may be determined element by element as 


1 n-l 
ie = : T 9: 


where IJ denotes a product over the range of 2 from the rows of the matrix AT 


Similarly a, is determined according to @p = apy - Ch_y Api Cy-4» Which 
is obtained from the rows of the matrix Aq}. 


6.0 COMPARISON OF METHODS 


Since the matrix which will be inverted is symmetric the preceding techniques will 
be somewhat simplified. 


6.1 Single Division 


It appears as if the single division scheme for S = D!D will yield the highest 
accuracy, since this method tends to make the leading diagonal of D closer to unity. 


For this method, the storage for n(n-1)/2 elements must be provided, since the 
determination of the elements of the Kth row of D-' utilizes all the elements of the 
rows £ =k which have already been determined. This disadvantage may be offset 
by suitable partitioning, but appears to become cumbersome. 


Furthermore, the required inverse A™ is not obtained until D-' is post-multiplied 
by its transpose. 


6.2 Partitioning 


This has the advantage that it enables a matrix of any size to be inverted, provided 
that suitable partitioning has been performed. However, it is apparent that this 
becomes rather cumbersome, and is not recommended unless a large number of zero 
or unit partitions exist. 


6.3 Bordering 


It has been shown that the bordering technique for a Schematic matrix may be 
obtained by the successive evaluation of algebraic expressions. These yield the 
elements of the inverse row by row, and computations are performed using rows 
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of the previously determined inverse one at a time. Provided that these rows are out- 
put onto tape (and read back as required) a maximum of (3n + 2) storage locations is 
required for the computations. 


However, the method is more time consuming than the single division scheme. 


6.4 Conclusion 


It appears that the bordering technique is the most appropriate method to use, 
on account of the simpler programming and storage requirements. 


The determination of A A~! and correcting of A7! by Hotelling’s method will 
probably have to be accomplished through a partitioning procedure. However, it is 
felt that it should be necessary to make only one iteration to obtain A™! with 
sufficient accuracy. 


In the event that a large number of iterations is necessary, the convergence 
may be improved by the methods described in Faddeeva Ch. 30. 
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